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ABSTRACT 

The Advanced Camera for Surveys on board the Hubble Space Telescope has been used to obtain deep 
high-resolution images of the giant early-type galaxy NGC 1316 which is an obvious merger remnant. These 
observations supersede previous, shallower observations which revealed the presence of a population of metal- 
rich globular clusters of intermediate age 3 Gyr). We detect a total of 1496 cluster candidates, almost 4 times 
as many as from the previous WFPC2 images. We confirm the bimodality of the color distribution of clusters, 
even in V — I, with peak colors 0.93 and 1.06. The large number of detected clusters allows us to evaluate the 
globular cluster luminosity functions as a function of galactocentric radius. We find that the luminosity function 
of the inner 50% of the intermediate-age, metal-rich ('red') population of clusters differs markedly from that 
of the outer 50%. In particular, the luminosity function of the inner 50% of the red clusters shows a clear 
flattening consistent with a turnover that is about 1.0 mag fainter than the turnover of the blue clusters. This 
constitutes the first direct evidence that metal-rich cluster populations formed during major mergers of gas-rich 
galaxies can evolve dynamically (through disruption processes) into the red, metal-rich cluster populations that 
are ubiquitous in 'normal' giant ellipticals. 

Subject headings: galaxies: star clusters — galaxies: elliptical and lenticular — galaxies: individual 
(NGC 1316) — galaxies: interactions 



1. INTRODUCTION 

Recent deep imaging studies of 'normal' giant ellipti- 
cal galaxies with the Hubble Space Telescope (HST) and 
large, ground-based telescopes have shown that these galax- 
ies usually contain rich globular cluster (GC) systems with 
bimodal colo r distributions (e.g., iKund u & Whitmore 2001; 
iLarsen et alJ I2OOI1) . Typically, roughly half of the GCs 
are blue, and half red. Follow-up spectroscopy with 8-m 
class telescopes revealed that both 'blue' and 'red' GC sub- 
populations are typically old ( > 8 Gyr, iForbes et alJ 1200 It 
ICphen. Blakeslee. & C6tdl l2003i) implving fliat the bimodal- 
ity is mainly due to differences in metallicity. The col- 
ors of the 'blue' GCs are usually similar to those of metal- 
poor halo GCs in the Milky Way and M31, while the mean 
colors of the 'red' GCs are sim ilar to those of the diffuse 
light of their host galaxies (e.^..lGeisler. Lee. & Kim 1996; 

^93- 



IForbes. Brodie. &' GrillmaiilfT997h . Hence, the nature of the 
'red', metal-rich GCs is likely to hold important clues to the 
star formation history of their host galaxies. 

One environment known to produce metal-rich GCs and 
bimodal color distributions is that of vigorous starbursts in- 
duced by mergers of gas-rich galaxies. Massive young GCs 
have been commonly found in merger s and young merger 
remnants using HST observations (e.g., Holtzman et al. 1992; 
ISchweizer 20_02.. and references therein). Follow-up spec- 
troscopy has confirmed the ages (a nd in one case even the 
high masses, i Maraston et al J 120041) of these young clusters 
predicted from their colors and luminosities (e.g.. Zepf et al. 
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ll995tLSchweizer & Seitzedll998h . Their metallicities tend to 
be near solar, as expected for clusters formed out of enriched 
gas in spiral disks. A natural interpretation of these data is that 
the metal-rich GCs in 'normal' giant ellipticals formed in gas- 
rich mergers at z > 1, and that the formation process of gi- 
ant ellipticals with significant populations of metal-rich GCs 
was si milar to that in galaxy mergers observed today (e.g., 
Schwe izeJ 1 987t IXshm an & Zepf 1992). 

However, one important, hitherto unsurmounted hurdle for 
this 'formation by merger' scenario has been the marked dif- 
ference in the luminosity functions (hereafter LFs) of old vs. 
young GC systems (e.g., Ivan den Bergh. 1995 ). The LFs of 
old GC systems of 'normal' galaxies are well fit by Gaussians 
peaking at A/,/ ^ —7.2 mag with a dispersion of cr ~ 1.3 mag 
(e.g., Whitmore Il997h . while young GC systems in merg- 
ers and young remnants have power-law LFs with indices of 
a ~ — 2 (e. g.. . Whitmore et al. 199 9). Recent calculations of 
dynamical evolution of GCs (including two-body relaxation, 
tidal shocking, and stellar mass loss) show that the least mas- 
sive GCs disrupt first as galaxies age, which can gradually 
transform the initial power-law LFs into LFs with Gaussian- 
like peaks or turnovers (e.g.. Fall & Rees 1977; BaumgardJ 
[^98; Fan & Zhang 2001; but see Vesperini 2001). Obser- 
vational evidence of this effect has been reported for a GC 
system near the center of M 82, featuring a very short GC dis- 
ruption timescale (de Grijs, Bastian, & Lamers l2003l) . 

Intermediate-age merger remnants with ages of 1-5 Gyr 
are ideal probes for studying the long-term dynamical effects 
on GC systems formed during a major merger. Such galaxies 
are still identifiable as merg er remnants through the ir morpho- 
logical fine structure (e.g., Schweizer & Seitzer 1992), yet are 
old enough to ensure that substantial dynamical evolution of 
the globular clusters has already occurred. 
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Recent HST studies of candidate intermediate-age merger 
remnants have revealed that their 'red' GC subpopula- 
tions show LFs consistent with power laws (as expected if 
formed during a rece nt merger event; Goudfrooij et al. 2001^ 
IWhitmore et al.1 '2002). However, the completeness limits of 
those studies were not faint enough to allow a detection of a 
turnover in the LFs. This Letter reports on new observations 
of the intermediate-age merger remnant NGC 1316 using the 
Advanced Camera for Surveys (ACS), installed on HST in 
March 2002, whose unprecedented sensitivity allows us, for 
the first time, to detect these turnovers. 

2. NGC 1316 

NGC 1316 is a prime example of an intermediate- age 
merger remnant. Extensive optical observations showed that 
NGC 1316 is a typical Morgan D-type galaxy with a surface 
brigh tness profile following an r"^/^ law (Schweizerl 119801 
Il981h . Its outer envelope includes several non-concentric 
arcs, tails and loops that are most likely remnants of tidal per- 
turbations, while the inner part of the spheroid is characterized 
by a surprisingly high central surface brightness and small ef- 
fective radius for the galaxy's luminosity. All of these features 
are consistent with NGC 1316 being the product of a dissipa- 
tive merger with incomplete d ynamical relaxation. Recently, 
Goudfrooij et al. (2001a; 2001b) discovered a significant pop- 
ulation of 3 Gyr old GCs of near-solar metallicity through a 
comparison of BVIJHK colors as well as Ha and Ca ll-triplet 
line strengths with population synthesis model predictions. 

In the follo wing we adopt a distan ce of 22.9 Mpc for 
NGC 1316 (see lGoudfrooii et al.l200ibl hereafter Paper I). 

3. OBSERVATIONS AND DATA ANALYSIS 

NGC 13 16 was observed with HST on March 4 and 7, 2003, 
using the wide-field channel of ACS and the F435W, F555W, 
and F8 14W filters, with total exposure times of 1 860 s, 6980 s, 
and 4860 s, respectively. In this Letter we concentrate on the 
F555W and F814W images, which reach significantly deeper 
than the F435W image. The final F555W and F814W images 
were constructed from long (18-21 min) exposures at 4 or 
6 dither positions, supplemented by a few short exposures to 
avoid saturation of the central regions. The individual images 
in each band were combined using the task MULTIDRIZZLE 
within IRAF/STSDAS v3.2. Saturated pixels in the long ex- 
posures were replaced by those in the short exposures while 
running MULTIDRIZZLE by setting the appropriate data qual- 
ity flag in the affected pixels of the long exposures. The high 
sensitivity of the ACS images allowed us to reach about 1.9 
mag deeper than the WFPC2 images used in Paper I did. 

WFPC2 images were taken in parallel to the ACS observa- 
tions, ^5'.5 away from the ACS pointings, in order to pro- 
vide estimates for the number density of compact background 
galaxies and foreground stars as a function of brightness and 
color (see Sect.|4}. The WFPC2 exposures consisted of a total 
of 5100 s in F555W and 3000 s in F814W. Each set of WFPC2 
images was also combined using MULTIDRIZZLE. 

The cluster selection procedure followed the one described 
in Paper I, to which we refer the reader interested in de- 
tails. Briefly, we first built smooth galaxy model images from 
the combination of an isophotal model composed of pure el- 
lipses and a median-filtered version of the image. The objects 
were then selected b y applyi ng the DAOFIND task from the 
DAOPHOT package ( Stetson 1987) to an image prepared by di- 
viding the F555 W image by the square root of the correspond- 
ing model image (thus having uniform shot noise characteris- 
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Fig. 1 . — Completeness functions used for the ACSAVFC photometry of 
NGC 1316 GC candidates. Lines (solid for V band, dashed for / band) are 
shown for five background levels: From left to right: 1600, 800, 320, 200, 
and 140 e^ pixel for an exposure time of 1100 s. 



tics over the whole image). The detection threshold was set 
at 3 cr above the background. Detections located within or 
on the edge of the dust features in NGC 1316 were excluded 
from the sample. As the resolution of the drizzled ACS im- 
ages is high enough to resolve typical GCs at the distance 
of NGC 1316, a further subselection was made by restricting 
a 'concentration index' (the difference in magnitudes mea- 
sured using radii of 2 and 5 pixels) to stay between 0.35 and 
0.75 mag. This effectively filtered out extended background 
galaxies as well as bright foreground stars and any remain- 
ing hot pixels. Aperture corrections were —0.49 (±0.01) mag 
in V and —0.56 (±0.01) mag in / for an aperture radius of 
3 WFC pixels (i.e., 0'.'15). These values were determined 
through measurements of several isolated, well-exposed GC 
candidates located throughout the field. The transformation 
of the instrumental magnitudes to Johnson/Cousins V and / 
was performed using the calibrations derived by the ACS In- 
strument Definition Team (Sirianni et al., in preparation). 

Completeness corrections were performed for five back- 
ground levels and several magnitude intervals by adding arti- 
ficial objects (in batches of 300) with a radial intensity profile 
derived from a fit to real GC candidates in each frame (using 
DAOPHOT routines). The object colors were set equal to the 
median color of GC candidates. Figure^shows the resulting 
completeness functions. 

4. TWO SUBPOPULATIONS AND THEIR LUMINOSITY FUNCTIONS 

The top panel of Figure|2depicts the V versus V~I color- 
magnitude diagram (hereafter CMD) for GCs in NGC 1316. 
The V — I colors of the GCs that are more luminous than any 
GC in the Milky Way are quite uniform and lie somewhat 
redward of the range covered by the MiUcy-Way halo GCs. 
Goudfrooij et al. (2001a) analyzed spectra of three of these 
bright GCs, and found all three to be ^ coeval with an age 
of 3.0 ± 0.5 Gyr and a solar metallicity (to within ±0.15 
dex). The fainter part of the CMD reveals a population of GCs 
with colors and luminosities that are consistent with t hose of 
halo GCs in our Galaxy (taken from the database of iHarrisI 
(199^, as well as a larger population with a mean color that is 
consistent with those of the aforementioned brightest, metal- 
rich GCs throughout the sampled range of luminosities. As 
to the potential effect of dust extinction within NGC 1316, 
our selection criteria seem to have rendered the number of 
GC candidates with colors redder than expected (given the 
photometry errors) negligible regarding the discussion below. 

The V — I color distribution of the GC candidates is shown 
in the bottom panel of Figure |2] for two brightness cuts. The 
bimodality appears more pronounced than in Paper I, and is 
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Fig. 2. — Top panel: V vs. V — I color-magnitude diagram for compact 
sources in NGC 1316. Circles represent GC candidates, while dots represent 
foreground star candidates, i.e. objects with concentration parameters (see 
text) that are consistent with that of the ACSAVFC point spread function. The 
inner 50% of the GC candidates are shown as open circles, while the outer 
50% are shown as filled circles. The hatched region represents the magnitude 
and color range for GCs in the Milky Way halo, placed at the distance of 
NGC 1316. Representative error bars are shown on the right-hand side of 
the diagram. The tilted line shows the reddening vector for Ay = 1 mag 
of extinction. Bottom panel: V — I color distribution for GC candidates in 
NGC 1316, showing all GCs with V < 2Q (open histogram), and GCs with 
V < 24.5 (filled histogram). The dashed lines represent GCs with V < 26 
before coiTection for contamination by compact background galaxies. 



clearly present down to = 26 where the typical uncertainty 
in V—I is ± 0. 12 mag. Prior to creating this diagram, the color 
histogram of candidate GCs was corrected for contamination 
by compact background objects. This correction was derived 
from the list of compact objects detected in the three WF chips 
of the parallel, drizzled WFPC2 images, selected according to 
the same criteria as the GC candidates in the ACS images. The 
color distribution of the objects in the WFPC2 images was 
smoothed with a three-bin kernel to diminish small-number 
noise and scaled to the area of the ACS image prior to sub- 
traction from the ACS color histogram. The effect of this cor- 
rection is shown in Figure|2] 

The luminosity functions of the blue and red GC subpop- 
ulations are shown in Figure |3] and were determined as fol- 
lows. We first separated the GC candidates into 'blue' (0.55 < 
V-I < 0.97) and 'red' (1.03 < V-I < 1.40) subsamples. 
We avoided the overlapping region (0.97 < V — I < 1.03) 
to minimize the mixing of the two subsamples. A complete- 
ness correction was then applied to each GC by dividing its 
count of 1 by the completeness fraction corresponding to its 
brightness and background level. The latter fraction was cal- 
culated from the functions shown in Fig.^ using bilinear in- 
terpolation in log(background) space. The LFs were further 



corrected for contamination by compact background objects 
using the smoothed and scaled LF of the objects in the par- 
allel WFPC2 images, as described in the previous paragraph. 
The corrected LFs are shown as solid lines in Fig.|3] whereas 
the raw, 'observed' LFs are shown as dotted lines. 

As found in Paper 1, the LFs of the (full) blue and red GC 
subsystems in NGC 1 3 1 6 are dramatically different from each 
other The shape of the LF of the blue GC system is consis- 
tent with that for GCs in 'normal' galaxies and peaks near 
the expected turnover magnitude of My = —7.2. In contrast, 
the LF for the red GC system is consistent with a power law 
{(j){L)dL oc dL with a = -1.75 ± 0.07) down to the 50% 
completeness limit of the data. At face value, these results do 
not provide significant evidence to suggest an evolution of the 
LF from a power law to a Gaussian. 

However, the disruption timescales of dynamical effects 
thought to be responsible for the Gaussian shape of the GC 
LFs of 'normal' galaxies through preferential depletion of 
low-mass GCs (e.g.. Fall & Zhang 2001, hereafter FZOl) de- 
pend on galactocentric distance. The stronger tidal field in the 
central regions does not only yield stronger tidal shocks than 
in the outskirts, it also imposes a more stringent zero-energy 
surface on the GCs so that stars can be more easily removed 
from the GCs through two-body relaxation. Hence, the ap- 
pearance of a turnover in the GC LF can be expected to first 
occur in the inner parts of galaxies (see also Gnedin 1997). 

With the 1496 GC candidates in our ACS images of 
NGC 1316, we are now in a position to test this idea for a 
^ 3 Gyr-old merger remnant with adequate statistical signif- 
icance. To this end, we divide the red GC system into two 
equal-size parts sorted by projected galactocentric radius. For 
our sample, this division ends up at 85'.'2, or 9.4 kpc from the 
galaxy center. The result is shown in Fig.|3lc) and (d): The 
LF of the outer half of the red GC system still rises all the way 
to the detection limit, but the LF of the inner half of the red 
GC system shows a clear flattening consistent with a turnover 
at ~ 1.0 mag fainter than that of the blue GC system. 

We compare these results with predictions of GC disrup- 
tion models by FZOl in Fig. |3je) and (f). The FZOl mod- 
els, whose initial GC velocity distribution involves a radial 
anisotropy which is similar to that of halo stars in the solar 
neighborhood, show that an initial power-law or Schechter 
mass function evolves in 12 Gyr to a peaked mass distribu- 
tion similar to that observed for the Milky Way GC system. 
To predict the shapes of luminosity functions as a function of 
age, we use the FZOl model that uses a Schechter initial mass 
function in combination with the stellar evolution models of 
■Maraston et aL (.2001) . We use Z/Zq = 1.0 for populations 
of ages 1.5, 3, and 6 Gyr, while we use Z/Zq = 0.02 for 
the 12 Gyr old population, similar to the median metallicity 
of halo GCs in our Galaxy. Prior to plotting, all resulting 
LFs were shifted by -1-0.2 mag so as to let the peak luminos- 
ity for the 12 Gyr old, metal-poor model be at Aly = —7.2. 
As Fig. 13 f) shows, the location of the turnover found for the 
inner half of the red GC system is in remarkably good agree- 
ment with those of the FZOl models for intermediate ages and 
solar metallicity. It thus seems reasonable to assume that the 
LF of this second-generation GC system will evolve further to 
become consistent with the LFs of red GCs in 'normal' giant 
ellipticals. 

5. CONCLUSION 

We have presented the discovery of a clear flattening con- 
sistent with a turnover in the LF of the inner half of the 
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Fig. 3. — y-band luminosity functions (LFs) of GC candidates in NGC 1316. Panel (a): LP of the full blue subsystem. Panel (b): LP of the full red 
subsystem. Panel (c): LP of the outer 50% of the red subsystem. Panel (d): LP of the inner 50% of the red subsystem. Dashed histograms [only in panels (c) 
and (d)] mark uncorrected (observed) LPs, dotted histograms mark LPs corrected for incompleteness, and solid histograms mark LPs corrected for completeness 
and background galaxies. Histograms are filled for magnitude bins brighter than the overall 50% completeness limit of the sample plotted, and open beyond 
it. The smooth curve in panel (b) is a power-law fit to the LP. The dotted vertical lines represent the predicted turnover magnitude for 'old' GC systems (i.e.. 
My = —7.2). Note the appearance of a flattening consistent with a turnover in the LP of the inner 50% of the red subsystem, which is not present in the outer 
50%. Panel (e): Same LP as in panel (a), but at a different scale and with formal errorbars. Panel (f): Same LP as in panel (d), but at a different scale and with 
formal errorbars. The smooth curves in panels (e) and (f) represent predicted LPs based on a combination of the Pall & Zhang (2001) GC disruption models 
and the population synthesis models of Maraston et al. ( 2001). Curves are drawn for the following populations: ages of 1.5 (short-dashed Hne), 3.0 (long-dashed 
line), 6.0 Gyr (dash-dotted line), and 12 Gyr (solid line in panel (e)). Z/Zq = 0.02 is assumed for the 12 Gyr curve, whereas Z/Zq = 1.0 is assumed for the 
three younger ages. The peak amplitudes of all model curves in panels (e) and (t) are scaled to coincide with those of the underlying histograms. 



red, metal-rich GC subsystem of NGC 1316, a ~3 Gyr-old 
merger remnant. The tur nover magnitude is consistent with 
predictions of thel Fall & Zhang ( ,2001) models at that age and 
metallicity. The implication is that the dynamical evolution of 
metal-rich GC populations formed in gas-rich galaxy merg- 
ers can change their properties so that they become consistent 
with the properties of the red, metal-rich GC subsystems that 
are ubiquitous in 'normal' giant ellipticals. This discovery 
weakens the one remaining argument against the 'formation 
by merger' scenario for the metal-rich GC subsystems in 'nor- 
mal', giant ellipticals, namely that the LFs of GC systems in 
young and intermediate-age merger remnants are inconsistent 



with the Gaussian LFs of old systems. Our results support a 
picture in which the formation process of giant ellipticals with 
significant populations of metal-rich GCs was similar to that 
in gas-rich galaxy mergers observed today. 
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